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Abstract 

Owing to the constraints on the fabrication of y-ray coding plates with many pixels, 
few studies have been carried out on y-ray computational ghost imaging. Thus, the 
development of coding plates with fewer pixels is essential to achieve y-ray 
computational ghost imaging. Based on the regional similarity between Hadamard 
subcoding plates, this study presents an optimization method to reduce the number of 
pixels of Hadamard coding plates. First, a moving distance matrix was obtained to 
describe the regional similarity quantitatively. Second, based on the matrix, we used 
two ant colony optimization arrangement algorithms to maximize the reuse of pixels in 


the regional similarity area and obtain new compressed coding plates. With full 


sampling, these two algorithms improved the pixel utilization of the coding plate, and 
the compression ratio values were 54.2% and 58.9%, respectively. In addition, three 
undersampled sequences (the Harr, Russian dolls, and cake-cutting sequences) with 
different sampling rates were tested and discussed. With different sampling rates, our 
method reduced the number of pixels of all three sequences, especially for the Russian 
dolls and cake-cutting sequences. Therefore, our method can reduce the number of 
pixels, manufacturing cost, and difficulty of the coding plate, which is beneficial for 
the implementation and application of y-ray computational ghost imaging. 

Key words: y-ray computational ghost imaging; Regional similarity; Hadamard coding 
plate; 

1. Introduction 

Computational ghost imaging (CGI), which is less expensive than traditional y-ray 
imaging, can reduce the scientific costs.'" By replacing the traditional detectors with 
bucket detectors and spatial modulators, CGI has attracted much attention as a new type 
of imaging.'*! Moreover, the crosstalk noise between the detectors and the radiation 
intensity significantly affects the image quality. CGI can effectively avoid crosstalk 
noise between detectors and reduce the radiation intensity required by imaging. In 
particular, it avoids the limitation of the size of detectors on image spatial resolution, 
which offers a new possibility for more precise radiation imaging. 

However, designing and manufacturing a suitable spatial modulator has always 
been a challenge for achieving CGI. The earliest breakthrough is visible light. A spatial 
light modulator (SLM) or digital micromirror (DMD) "l has been used to modulate 
the light field. However, neither the SLM nor the DMD can modulate ray particles with 
certain penetration ability. In the field of rays, the only way to achieve this is to take 
advantage of the attenuation of rays by coding plates. 

The difficulty of fabricating coding plates significantly delays the realization of 
X-ray and neutron ghost imaging. In 2020, He et al.'© 7! used electroplating to fabricate 
Hadamard coding plates made of Au and photochemical etching to fabricate Hadamard 
coding plates made of Cu. Their pixel sizes were 10 x 10 and 150 x 150 um, and their 


thicknesses were 10 and 150 um, respectively. Using the two coding plates as spatial 


modulators, they achieved X-ray CGI with resolutions of 10 and 150 um. In a 
subsequent study of neutron CGI, He et al.'®! used the Bosch process for deep silicon 
etching and Gd2Oz3 powder filling to produce a thermal neutron coding plate with a 
thickness of 300 um. They also conducted the first thermal neutron CGI experiment on 
the Dongguan spallation neutron source device. 

However, few studies on y-ray CGI are available at this moment. This is because 
the penetration ability of y-rays is generally higher than that of X-rays and thermal 
neutrons. For example, if the energy of y-rays is above 100 keV, the thickness required 
for high Z-encoding plates must be more than 1 mm, which is currently impossible to 
achieve using etching processes. Similar devices are also needed in y cameras to correct 
y rays, such as those in the uniformly redundant arrays coding plate (URA) and 
modified uniformly redundant arrays coding plate (MURA). '*! These coding plates are 
made by cutting the metal into N small blocks. Each small metal block is then placed 
into a coding plate according to the numerical distribution of the respective correction 
matrices. However, to achieve the same imaging resolution, N° pixels are required in 
the CGI coding plate. It is difficult to employ this method to create coding plates with 
many pixels. 

In addition, to achieve a higher imaging resolution, the number of required 
measurements has increased dramatically, and the number of subcoding plates required 
has also increased significantly, resulting in a very large number of pixels on the coding 
plate. This eliminates the cost advantage of a single-pixel imaging system, and the cost 
savings from detectors are wasted on coding plate manufacturing. For a high-quality y- 
ray imaging system, it is difficult to create coding plates with a large number of pixels. 
Therefore, the development of coding plates with fewer pixels is essential to achieve y- 
ray CGI. 

To date, few good solutions to this problem have been reported. The reasons for 
this are as follows: (1) relevant research is insufficient because of the short history of 
ray CGI; (2) in the optical field, it is easy to change the distribution of modulated light 
fields because of the widespread use of SLM and DMD. The study of light-field 


modulation mainly focuses on achieving high-quality imaging under sampling. 


If the undersampling method is used in y-ray CGI, the number of required coding 
plate pixels can also be reduced. Nevertheless, this is still excessive for y-ray CGI. It 
will be easier to fabricate coding plates that are suitable for y-rays if the number of 
pixels is reduced further. 

Compressed sensing (CS) theory provides an excellent theoretical framework for 
image reconstruction in CGI with undersampling. It often uses some special measured 
patterns, such as the Hadamard mask and random mask, and iterative optimization 
methods, such as the L1-magic algorithm, orthogonal matching pursuit (OMP) 
algorithm, and TVAL3 algorithm, for image reconstruction.'°! The undersampling 
method can reduce the number of pixels required by the coding plates, but this also 
means a higher standard in the design of coding plates.'!""°! That is to say, in the limited 
number of measurements, the coding plate should maximize the degree of "fluctuation" 
of the measurement data so that valid information can be retained as much as 
possible.''*?°! Random masks have the advantage of being easy to manufacture. One 
effective method to reduce the sampling rate further is to use orthogonal basis types of 
masks, such as Hadamard coding plates, to remove the redundancy in random 
measurements. 

Currently, Hadamard coding plates are widely used in CGI. Hadamard coding 
plates are adopted because the Hadamard matrix has some interesting features, such as 
good orthogonality and the ability to remove redundant information in natural scenes.'"7! 
Then, the question arises, which subcoding plates should be selected to form a new 
plate as small as possible and able to maintain high-quality imaging? Researchers have 
proposed various undersampling measurement sequences based on the Hadamard 
matrix. They hope that the careful design of multiplexing matrices will result in 
potential performance improvements in CGI. 

Sun et al.'°! proposed a Russian doll measurement sequence in 2017. This 
sequence was designed according to both the characteristics of the numerical 
distributions and the total number of subcoding plate blocks. This sequence can achieve 
good recovery results, with a sampling rate of only 16%. In 2019, Yu et al."!®! found 


that rearrangement of the measurement sequence of Hadamard coding plates might 


change the coherent area of each imaging point, which might in turn lead to a change 
in the quality of imaging. If the measurement sequence of the coding plates is arranged 
in ascending order in terms of the number of blocks, the coherent region of each 
imaging point will naturally become smaller, resulting in a high-quality imaging 
recovery at smaller sampling rates. This measurement sequence, known as the "cake- 
cutting" sequence, has a better effect than the Russian dolls and random Hadamard 
measurement sequences in undersampling. In addition, because the image 
reconstruction effect of the compressed sensing algorithm is closely related to the 
selected sparse matrix, some researchers have proposed measurement sequences based 
on different sparse matrices, such as the Harr sequence.''®: 7°! At a sampling rate of 25%, 
this Harr sequence can better restore the image characteristics. 

In 2020, Xiao proposed two new measurement sequences from the perspective of 
graph decomposition: the total variation (TV) method and the total wavelet transformed 
coefficients (TW) method.'*"! The subcoding plates are reordered through their TV and 
TW in ascending order to achieve the best performance. The numerical simulation and 
experimental results show that the TV, CC, and TW orders can provide almost the same 
reconstruction at deep compressive sampling, and even better performance at a 
sampling ratio above 30%. 

Yu proposed a concept called selection history, which can record the Hadamard 
spatial folding process, and built a model based on it to reveal the formation 
mechanisms of different orderings and deduce the mutual conversion relationship 
among them.'"”! Based on this, a weight ordering of the Hadamard basis was proposed. 
Both the numerical simulation and experimental results have shown better 
reconstruction quality with a lower sampling rate than traditional sorting methods. 


Many orderings have emerged, but their relations remain unclear and lack a unified 


theory to explain the inherent math and physical mechanisms. In particular, there is a 
each undersampling sequence. In addition, most researchers have concentrated on how 


to compose a multiplexing matrix in the undersampling method to reduce the number 


of pixels of the coding plates (mask). However, the method that directly reduces the 


number of pixels of the coding plates by changing the connection between the 

subcoding plates is ignored. 

In this study, on the premise that the size of coding boards cannot be reduced on a 
large scale owing to the unification of interpretation theories at present, the internal 
correlation of the subcoding plate is studied from another perspective, and the new 
connection mode between subcoding boards is studied to reduce the number of pixels 
of coding plates. 

The specific approach is as follows: we propose an optimization algorithm that 
can create “compressed coding plates.” A compressed coding plate is generated by 
compressing a traditional coding plate. The algorithms satisfy the requirements of both 
full sampling and undersampling. This study consisted of three parts. 

1) The characteristics of Hadamard coding plates were studied, and a moving distance 
matrix was obtained to quantify the regional similarity between the subcoding 
plates. 

2) Two ant colony optimization arrangement algorithms were designed to reuse pixels 
of the same region to the largest degree between the subcoding plates. 

3) The optimum method was tested, and its outcomes were discussed in terms of full 
sampling and undersampling. 

2. Methods 

2.1 Computational ghost imaging 


The schematic of a common y-ray CGI device is shown in Fig. 1. 


Fig. 1 Conceptual drawing of the experimental setup used for y-ray CGI. The left and right 
parts of the image show the overall layout and a partial enlargement of the equipment, 


respectively. The collimator (1) ensures that the particles produced by the y-source are only 


incident on the imaging area; the coding plate (mask) (2) is composed of all the subcoding plates, 
whose pattern is encoded on metal; the particles are modulated on the subcoding plate (3), and the 
object can be placed between the collimator and subcoding plate; the signal-pixel detector (4) 
collects the intensity of rays; a computer (5) is used to reconstruct the images. 

In CGI, it is necessary to move the coding plate after each measurement to change 
the shape of the mask in the imaging area and thus, the distribution of the ray field. In 
this study, the mask of the corresponding part of each measurement is called a 
subcoding plate. The measurement process can be expressed using Eq. 1: 

i= Pi*O (1) 
where 

O represents N x N pixel values in an object and is transformed into a column 
vector. 

Pi represents N x N modulation values corresponding to the itn subcoding plate (a 
measured pattern) and is transformed into a row vector. 


Ii represents a measured value obtained from the in measurement. 


In CGI, the measured value J, is mathematically equivalent to the inner product 


of each pattern Pi and object image O. The spatial information of objects is acquired by 
a detector and subcoding plates, not by array detectors. 


Traditional ghost imaging is a well-known imaging technique. Its definitions are 


as follows. An object image O, can be reconstructed using P and I,. 
O, *(x,y) = (LP (x, yY)}- (Ui) (By) (2) 
where (.) denotes the average of all the CGI measurements. 


With the development of ghost imaging technology, some researchers have 
reconstructed images from the perspective of solving equations. For example, Zhang 


21 and Xue proposed singular value 


proposed pseudo-inverse ghost imaging 
decomposition ghost imaging.'*?! Then, optimization algorithms based on compressed 


sensing (CS) have been employed to reconstruct images under undersampling, such as 


the orthogonal matching pursuit (OMP) and TVAL3 algorithms. 

In addition, some patterns, such as the Gaussian, Bernoulli, DCT, and Hadamard 
matrices, have been used in CGI. The Hadamard matrix can generate better-quality 
images compared to other measurement matrices.'!7! Therefore, the Hadamard coding 
plate is widely used in CGI. 

2.2 Regional similarity 

Because of the uncorrelated character of any two columns or rows, the Hadamard 
matrix, as a modulation orthogonal matrix, can effectively avoid correlation noise 
among the pixels, and the reconstructed image quality can be greatly improved at the 
same sampling rate.'!*! The Hadamard matrix is generally constructed recursively from 
low to high order, as shown below: 


H= Hoen Hyen 
2" H H (3) 


y(n)» = 20) 


The traditional Hadamard coding plate comprises K x K-block subcoding plates. 
The molding method of the subcoding plates converts the 1D code elements into 2D 
patterns, as shown in Fig. 2(a—c). 

Owing to the recursive property of the Hadamard matrix, it can be observed that 
the front part of a row in the matrix is the same as the back part of another row (Fig. 
2(a)). This rule causes the left area (Fig. 2(b)) of one subcoding plate to be completely 
consistent (Fig. 2(d)) with the right area (Fig. 2(c)) of another subcoding plate. This 


phenomenon is known as regional similarity between two subcoding plates. 
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Fig. 2 Modeling method of subcoding plates and regional similarity between subcoding 
plates. (a) Partial numerical distributions of the Hadamard matrix with a size of 64° order; (b) and 
(c) represent two subcoding plates forming the 10m and 42na rows of the matrix, respectively, 
which are clearly marked by purple and green borders, and the regions of same pixel values are 


marked in red; the two subcoding plates were compressed according to the regional similarity (d). 


Table 1 Moving distance matrix between all subcoding plates 


Note: The subcoding plate corresponding to the i-row of the Hadamard matrix is called the No.i 


subcoding plate C;. The value of regional similarity represents the number of columns in the same 
part between the left area of Cm and the right area of Cn. 


If we want to use regional similarity to compress the coding plate, it is necessary 
to calculate the degree of regional similarity between subcoding plates quantitatively. 
To achieve this goal, a moving distance matrix M is employed to describe the regional 
similarity. The moving distance matrix Mes corresponding to the Hadamard matrix 
with a size of 64? (64 subcoding plates) is presented in Table 1. 

According to the matrix in Table 1, the regional similarity between the subcoding 
plates has the following characteristics. 

1) Approximately 90% of the elements in the moving distance matrix take the value 
of 8, which indicates that there is no regional similarity between most subcoding 
plates. However, there are 320 pairs of subcoding plates with regional similarities, 
which indicates that the aim of reducing the number of pixels of the coding plates 
can be achieved by adjusting the order of the subcoding plates. 

2) In the matrix, there is a sporadic continuous oblique distribution of low-value 
elements, which indicates that there are interrelated characteristics between 
subcoding plates with regional similarities. This makes it possible to use 
continuous subcoding plates with regional similarities. 

Therefore, based on the moving distance matrix, we reuse the same region between 
the subcoding plates to reconstruct a new coding plate with fewer pixels, as shown in 
Fig. 2(d). The new coding plate can compress a traditional coding plate and is called a 
compressed coding plate. 

2.3 Ant colony optimization arrangement algorithms 

2.3.1 Ant colony algorithms 

To maximize the use of regional similarity between the subcoding plates and thus 
minimize the pixels of the compressed coding plate, it was necessary to redesign a 
measurement sequence of the subcoding plates according to the moving distance matrix 


obtained in the previous section. The problem is modeled as follows. 


N-1 
Min S=9 mii, j,) 
=1 


t 


J, = la (4) 
i=, 

St. i [1,2,3...,N] 
i #iand { j|je[1,2,...,1) O (ii, +1,.....N]} 


where N represents the number of points in a set; mij represents the distance from point 
i to point j, and the distance between points can be expressed by moving matrix M; S 
represents the sum of all distances in a path; t represents the tm point in a path, and it 
ensures that the distance between adjacent points is calculated. 

In the model, the goal is to find a path that minimizes S; the constraints are as 
follows: 1) each point can only pass once; and 2) all points need to be passed. This 
problem is similar to the travelling salesman problem (TSP). It is a mathematical 
combinatorial optimization problem and an N-P problem. The ant colony algorithm is 
a mature and effective algorithm for solving the TSP and other combinatorial 
optimization problems.'?*! In this study, each subcoding plate was used as the node of 
the ant colony path space, and the moving distance matrix was used as the node spacing. 
The path of an ant passing through all nodes is a feasible solution to the problem, and 
all the paths of the entire ant population constitute a feasible solution space. 

Affected by the communication mechanism of ants, their release of pheromones, 
and the volatilization mechanism of pheromones, more pheromones will be released on 
the short path (the optimal solution), whereas the pheromones on other solutions will 
volatilize slowly with iterations. Ants prefer paths with high pheromone concentrations. 
As the number of iterations increases, more pheromones remain on a shorter path and 
more ants choose this path. Eventually, under the action of positive feedback, ants 
choose the best path, and the corresponding subcoding plate sequence is the optimal 


solution to this problem. 
In addition, when ant k is at point t, the probability (Pi (t)) of choosing the next 


point j from point i mainly considers two factors: first, the distance between points i 


and j, which is called the heuristic factor 77,(t). Second, the feasibility from points i 


and j, which is called pheromone concentration r,,(t). P (t) can be calculated using 


the following formula: 
[aof no] 
Zao! two! 
PHs 0, j £ jy 


€ jy (5) 


where q represents the importance of pheromones, / represents the importance of the 
heuristic factor, and j,(t) is a selectable point set. 


After obtaining the path of the generation ant, the path is used to update the 


pheromones. The updated formula is as follows:. 
7, (t+) =(-p)z,(t)+ Az, 


Nant 
AT = LAT; (6) 
n=1 
AT; = a 


where p is the volatilization coefficient of pheromones, N „ is the number of 
generation ants, Az, is the sum of pheromones left by ants on the path between point i 
and point j, Ar, is the pheromone left by the kin ant on the path between point i and 


point j, Q, is the total pheromone of an ant, and L, is the total length of the path taken 


by the kun ant. 
2.3.2 Method for optimum coding plate 

According to the algorithm described in Section 2.3.1, the measurement sequence 
of the subcoding plates under the maximum use of regional similarity can be calculated. 
However, the coding plate manufactured according to the sequence is very long. 
Considering the actual production and usage, we need to optimize the coding plate 
further. A compressed coding plate should be designed, which should have a rectangular 


shape. Thus, it is called a rectangular compressed coding plate. This means that there 


should be no significant difference in terms of the length of each segment of the coding 
plate, and the total number of pixels of the compressed coding plate should be 
minimized. To solve the above problems, two methods can be employed for obtaining 
the optimum coding plate: the one-objective ant colony segmentation algorithm 
(OACSA) and the multi-objective ant colony reconstruction optimization algorithm 
(MACRA). 

(1) One-objective ant colony segmentation algorithm 

The core idea of this method is that the long measurement sequence obtained by 
the ant colony algorithm is divided into multiple short measurement sequences by 
selecting the appropriate segmentation points, and the long measurement sequence does 
not need to be changed. These short measurement sequences are reconstructed to obtain 
the corresponding rectangular compressed coding plates (supplemental pixels and 
parallel combinations). The reconstruction process may cause an increase in the total 
number of pixels in the coding plate. After the analysis, the position of the segmentation 
point should meet the following three criteria: 

1) The length of each segment is approximately equal. 

2) The segmentation points should be carefully selected to preserve the 
completed subcoding plate to avoid unnecessary pixel supplementation. 

3) The segmentation points should be selected between the subcoding plates with 
low regional similarity to ensure that the number of supplemental pixels involved in the 
reconstruction process is small. 

Based on these three criteria, this study proposes a segmentation-point selection 
method for a single-sequence coding plate, and its flowchart is shown in Fig. 3. 
According to the segmentation point selected by the algorithm, the single-sequence 
coding plate is divided to obtain multiple coding plate segments. Then, multiple 
segments are reconstructed to obtain the corresponding compressed coding plate 


(supplemental pixels and parallel combination). 


Settings of parameters M, L 
and J 


| 


q=1 4 


Allowable length range: 
[L-J, L+J] 
Criteria 1 l 


Nq=Nq-1+L 


4 


Optional segment point set: 
A=[Nq-J,Ng+J] 


Ss Optional segment point set: 
Criteria 2 B 
(Criteria 3 Ideal segment point C from B | 


$ 


Nq: the position of C | 


qt+ 


| M coding plate segments | 


Fig. 3 Segmentation-point selection method for a single-sequence coding plate 


Note: M represents the number of segments, L represents the mean length of each segment, J 
represents the allowable length deviation in each segment, Nq represents the location of the qin 


ideal segmentation point, and No = 0; M represents the number of segments and jy = floor( Lay : 
K 


Lan represents the length of a single-sequence compressed coding plate; K is the order of the 


subcoding plate. 


(2) Multi-objective ant colony reconstruction optimization algorithm 

The reconstruction of a rectangular compressed coding plate is a problem of 
reconstructing multiple independent segments and is an optimization problem with 
constraints. To solve this problem, we included the constraint (each segment is 
essentially of the same length) as the penalty part in the objective function, as shown in 


Eq. 7. 


N seg 


Min S=%s, +y* max (S, Sz, 59-445, ) 
k=1 se 


(7) 
Nx 

where : s, = > m(i, j) 
t=1 


where N, is the total number of segments and also represents the number of ants. Sp 


is the distance traveled by the kin ant in the same group, and it represents the length of 


the km segment, N, represents the number of points that the ant needs to cross, and y 


is the penalty factor coefficient. 

The MACRA was employed to solve the model. Through continuous iterative 
optimization, the total distance was gradually reduced, and the gaps between the lengths 
of all segments were bridged in the multi-objective colony algorithm. 

2.4 Evaluation index 

For a quantitative measurement of the degree of reduction in the number of pixels 
of the coding plate with ant colony optimization arrangement algorithms, this study 
defined a quantitative index-compression ratio, a. It represents the ratio of the number 
of pixels between the compressed and traditional coding plates. The smaller the value, 


the better the compression effect. Its formulas are as follows: 


rp 


a= i=1:M 


K’ 


p- 5 dK xk 


j=1:pi-1 


(8) 


where Di represents the number of pixels for the in segment, K represents the order of 
a subcoding plate, Pi represents the number of subcoding plates in each segment, di, 
denotes the minimum movement distance from the jm to the (j+1)th subcoding plate in 
the im segment, the value of i varies from 1 to M, and M is the number of segments. 
3. Results and discussions 

We applied the ant colony optimization arrangement algorithms described in 
Section 2 to reconstruct the full sampling traditional Hadamard coding plate and obtain 
a compressed coding plate. 


Subsequently, to demonstrate the feasibility of undersampling, we obtained the 
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compression ratios of compressed coding plates corresponding to the Harr, cake-cutting, 
and Russian dolls sequences at different sampling rates. The compression effects were 
compared and analyzed. 
3.1 Compressed coding plate in full sampling 

For all the subcoding plates corresponding to the 1024-order Hadamard matrix, 
the regional similarity between the subcoding plates was quantified using the moving 
distance modeling method to obtain the moving distance matrix. The OACSA and 
MACRA (the main parameters of the two algorithms are listed in Table 2) were used to 
reconstruct the rectangular Hadamard compressed coding plates, which are shown in 


Fig. 4. The compression results are shown in Fig. 5. 


(a) (b) 


Fig. 4 Traditional and rectangular compressed coding plate. 
(a) Traditional coding plate in full sampling, (b) compressed coding plate in full 
sampling. 


Table 2 Main parameters of the ant colony optimization arrangement algorithms 


Parameter Colony Iteration Importance of Pheromone Heuristic factor Segment 


name size number pheromones evaporation rate coefficients number 


Value 10240 10240 1 0.5 1 16 


1200000 - 35 p Minimum distance of movement 
B Number of pixels 
1048576 


The value is 32, accounting for 59.9% 


617611 


1000000 + 
Decrease oe 
800000 + 45.8% 41.1% 
600000 + 568328 
400000 + 
200000 + 
0 


Number of pixels 


Minimum distance of movement 


u 
Tf 


> 


1 mn fi 1 
200 400 600 800 1000 


Tradition OACSA MACRA 0 
Type of coding plate Sub-coding plate number 
(a) (b) 


Fig. 5 Experimental results in full sampling. (a) Compression results of the ant 
colony optimization algorithms in full sampling; (b) minimum movement distance 


between adjacent subcoding plates in the single-sequence coding plate. 


As shown in Fig. 5(a), both the OACSA and MACRA can significantly reduce the 
number of pixels, and their corresponding compression rates a are 54.2% and 58.9%, 
respectively. Moreover, the OACSA has a better compression effect. The reasons for 
this are as follows: 

With constraints, the MACRA is the basic solution to this optimization problem of 
combination, but it is more complicated than the OACSA. The MACRA requires a long 
time to obtain an answer. The answers are easily affected by the original figures. We 
believe that the OACSA may be a better solution to this challenging problem. The 
OACSA was applied in two phases: 

Phase 1: Without constraints (to make the coding plates rectangular), the one-way 
ant colony algorithm in the OACSA can solve the optimization problem. It is relatively 
simple and single-sequence coding plates with very few pixels can be made. 

Phase 2: When constraints are involved, single-sequence coding plates are 
segmented and combined to create rectangular compressed coding plates. If appropriate 
segmentation points are selected through the segmentation point selection method, the 
entire process results in only a minor increase in the number of pixels. 


The minimum distances of movement between adjacent subcoding plates in the 


single-sequence coding plate generated from the OACSA are shown in Fig. 5(b). It can 

be observed that the points with a maximum value of 32 and minimum moving distance 

account for 59.9%, and their distribution is wide. This phenomenon is favorable for the 

segmentation point selection method; thus, only a small number of pixels are added. 

This means that there is no significant difference between the number of pixels in the 

rectangular compressed coding plate and the single-segment arrangement scheme. Thus, 
we conclude that the OACSA is better than the MACRA. 


3.2 Compressed coding plate in undersampling 


In the previous section, it was proven that both optimum algorithms have good 
compression effects for the Hadamard coding plate in full sampling; however, the 
undersampling compression effects of the algorithm are unknown. Therefore, in this 
section, for the Harr, Russian dolls, and cake-cutting sequences, the OACSA, which 
has a better compression effect, is applied to reconstruct compressed coding plates at 


sampling rates of 5%—40%. The a values of the compressed coding plates are shown in 


Figure 6(a). 
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Fig. 6 Experimental results in undersampling. (a) Compression ratio of three 


sequences at different sampling rates; (b) trends of the number of blocks of the sub- 


coding plates in different sequences. 


As shown in Fig. 6(a), the ant colony optimization arrangement algorithm has a 
certain compression effect on the three sequences at different sampling rates, especially 
the Russian dolls and cake-cutting sequences. The reduction in the number of pixels in 
the two sequences can reach 40%. 

According to relevant literature!’® 181, the Russian dolls sequence can achieve a 
good imaging result at a sampling rate of 25%, that is, the number of pixels required 
for the coding plate is 25% of the number of pixels of the traditional Hadamard 
encoding plate. If the coding plate is compressed using the optimization arrangement 
algorithm, the number of pixels can be further reduced and the same imaging effect can 
be achieved without changing the modulation matrix. Therefore, the use of both the 
undersampling method and the ant colony optimization arrangement algorithm will 
significantly reduce the number of pixels in the coding plate. This will help lower the 
difficulty and cost of manufacturing the encoding plate and consequently build the 
foundation for the realization of y-ray CGI. 

In addition, the compression effect of the cake cutting sequence is the greatest, 
followed by the effect of the Russian dolls sequence, whereas that of the Harr sequence 
is not obvious. This result shows that the ant colony optimization arrangement 
algorithm has different compression effects on different types of sequences. The main 
reasons for this are as follows: 

A large area of the same pixel values is likely to appear between the subcoding 
plates with a small or similar number of blocks, that is, regional similarities can be 
easily found between such subcoding plates. 

The cake-cutting sequence is arranged in ascending order of the number of blocks. 


Russian dolls sequences are designed to take advantage of the fact that higher-order 


Hadamard matrices (H.,,,, ) contain low-order Hadamard matrices ( H m> m<n) . In 


their specific sorting method, patterns of subcoding plates with scaled versions of low- 
order matrices will be preferentially measured. Then, subcoding plates with a low 
number of blocks will also be preferentially measured. This method results in the cyclic 


growing trend of the number of blocks of Russian dolls. The Harr sequence is arranged 


in ascending order of the Harr value of the subcoding plate. According to these three 
different sorting rules, the trends of the number of blocks of the subcoding plates of 
different sequences are shown in Fig. 6(b). Both the cake-cutting and Russian dolls 
sequences are sorted in connection with the number of blocks, whereas the Harr 
sequence is not. As a result, the compression effects of the cake-cutting and Russian 


dolls sequences are much better than that of the Harr sequence. 


4. Conclusion 

In this study, optimum ant colony algorithms based on regional similarity are 
proposed to reduce the number of pixels and reconstruct a rectangular compressed 
coding plate. Without sacrificing the imaging quality, their compression rates can reach 
60% in full sampling. In undersampling, the proposed algorithm has a certain 
compression effect on each of the three sequences at different sampling rates, and the 
effects are obvious on the Russian dolls and cake-cutting sequences. The use of both 
undersampling and the proposed algorithm significantly reduces the number of pixels 
in the coding plate. This will help reduce the difficulty and cost of manufacturing the 
encoding plate and build the foundation for the realization of y-ray CGI. Finally, in 
addition to y-ray CGI, the proposed algorithms can also be applied to other types of 
CGI that require modulation by mechanical movement of the coding plate. 
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